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ABSTRACT 

We derive light curves of the afterglow emission from highly collimatcd jets if the power-law index 
(p) of the electron energy distribution is above 1 but below 2. We find (1) below the characteristic 
synchrotron frequency, the light curve index depends generally on p. (2) As long as the jet expansion is 
spherical, the light curve index above the characteristic frequency increases slowly as the spectral index 
of the emission increases. (3) Once the jet enters the spreading phase, the high-frequency emission flux 
decays as cx £~(p+ 6 )/ 4 rather than cx t~ p . All these results differ from those in the case of p > 2. We 
compare our analytical results with the observations on the GRB 010222 afterglow, and conclude that 
the jet model may be unable to explain the observed data. Thus, a more promising explanation for this 
afterglow seems to be the expansion of a relativistic fireball or a mildly collimated jet in a dense medium. 

Subject headings: gamma-rays: bursts — relativity — shock waves 



1. INTRODUCTION 

Gamma-ray burst (GRB) afterglows are believed to be 
emitted from a relativistic shock wave expanding in its 
surrounding medium via synchrotron radiation or inverse 
Compton scattering (ICS) of accelerated electrons in the 
shocked matter (Piran 1999; van Paradijs, Kouveliotou & 
Wijers 2000; Cheng & Lu 2001). To interpret the abun- 
dant data of afterglows, the effects of environments such as 
pre-burst stellar winds (Dai & Lu 1998; Meszaros, Rees & 
Wijers 1998; Chevalier & Li 1999, 2000) and dense media 
(Dai & Lu 1999, 2000; Wang, Dai & Lu 2000) have been 
discussed. On the other hand, jets are of particular inter- 
est because they have important implications on almost 
all aspects of the GRB phenomenon, e.g., the total energy 
that is released in an explosion, the event rate, the physi- 
cal ejection mechanism and the afterglow decay rate. The 
most exciting implication is that the transition of a rela- 
tivistic jet to the spreading phase can result in steepening 
of the afterglow light curve to the flux oc t~ p , as analyzed 
by Rhoads (1999) and Sari, Piran & Halpern (1999). Fol- 
lowing the analytical work, many numerical calculations 
have been performed and they are essentially consistent 
with the analytical results (e.g., Panaitescu & Meszaros 
1999; Moderski, Sikora & Bulik 2000; Kumar & Panaitescu 
2000; Panaitescu & Kumar 2001; Huang et al. 2000a, b, 
c; Wei & Lu 2000). The jet model seems to account well 
for a few well-observed afterglows with light curve breaks, 
e.g., GRB 990123 (Kulkarni et al. 1999; Castro-Tirado et 
al. 1999; Fruchter et al. 1999), GRB 990510 (Harrison et 
al. 1999; Stanek et al. 1999), GRB 991216 (Halpern et al. 
2000), GRB 000301C (Rhoads & Fruchter 2001; Masetti 
et al. 2000; Jensen et al. 2000; Berger et al. 2000; Sagar 
et al. 2000), GRB 000418 (Berger et al. 2001), and GRB 
000926 (Price et al. 2001; Harrison et al. 2001; Sagar et 
al. 2001a; Piro et al. 2001). 

GRB 010222 is the latest well-observed burst, whose 
optical-afterglow light curve has an earliest sharp break 
(Masetti et al. 2001; Stanek et al. 2001; Sagar et al. 



2001b; Cowsik et al. 2001). A popular explanation is that 
this afterglow might have come from a highly collimated 
jet with a flat-spectrum electron distribution (1 < p < 2). 
In this Letter we derive light curves of the emission when 
such jets expand in an interstellar medium (ISM) or in a 
stellar wind, and find that the jet model may be inconsis- 
tent with the afterglow data of GRB 010222. 

2. LIGHT CURVES 

Let's assume an adiabatic relativistic jet with an initial 
half opening angle of 0q, a laterally-spreading velocity of c s 
and a bulk Lorentz factor of 7. This assumption is valid if 
the energy density of the electrons accelerated by a shock, 
produced by the interaction of the jet with its surround- 
ing medium, is a small fraction e e of the total energy den- 
sity of the shocked medium or if most of the electrons are 
adiabatic, i.e., their radiative cooling timescale is larger 
than that of the jet expansion (Sari, Piran & Narayan 
1998). The energy density carried by magnetic fields is as- 
sumed to be another fraction es of the total energy density 
of the shocked medium, and thus the magnetic strength 
B = [327res7 2 'T-(r)TOpC 2 ] 1 / 2 , where m p is the proton mass 
and n(r) is the proton number density of the surrounding 
medium at shock radius r. We adopt a power-law density 
profile: n(r) = Ar~ s , where A = rt* x 1 cm~ 3 for the ISM 
(s = 0), and A = 3 x 10 35 A* cm -1 for the wind (s = 2) 
(Chevalier & Li 1999). 

We consider synchrotron radiation of the electrons ac- 
celerated by the shock. To calculate the spectrum and 
light curve, one needs to determine three break frequen- 
cies: the self- absorption frequency (is a ), the characteristic 
frequency {v m ) and the cooling frequency (v c )- The lat- 
ter two frequencies can be directly derived from the min- 
imum Lorentz factor 7 m and the cooling Lorentz factor 
7c which appear in the energy distribution of cooled elec- 
trons. As usual, we adopt a power-law injection of elec- 
trons with the energy distribution (just behind the shock 
front) given by dne/dje oc 7~ p for j m < j e < 7a/, where 
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7m = [3e/(^o'Ti?)] 1 ' 2 is the maximum electron Lorentz 
factor, which is calculated by assuming that the accelera- 
tion time equals the synchrotron cooling time. Here £ <~ 1 
is the ratio of the acceleration time to the gyration time, 
e is the electron charge, and <jt is the Thomson cross sec- 
tion. According to this electron energy distribution and 
the jump conditions for a relativistic shock, the electron 
number density and energy density of the shocked medium 
can be written as two integrals: J^ M (dn e /dj e )d-y e = Ajn 

and J 7M (-f e m e c 2 )(dn e /d'j e )d'y e — A"f 2 nm p c 2 e e , where m e 
is the electron mass. Such integrals combined with the 
assumption of a flat electron spectrum (1 < p < 2) lead to 

Equation (1) is different from the frequently- used mini- 
mum Lorentz factor -f m — [(p — 2)/(p — l)](m p /m e )e e j, 
which is also derived from these integrals for p > 2. Con- 
sequently, we will obtain expressions for afterglow light 
curves that differ from derived by Sari et al. (1998, 1999). 
In addition, the age of the jet could also provide a limit 
on -f m through constraining 7m in equation (1). However, 
we have found that the value of 7m inferred from the age 
limit is usually much larger than the one from the limit 
that the acceleration time equals the synchrotron cooling 
time, and thus the age limit on j m can be ignored. The 
7c the Lorentz factor of electrons that cool on the ex- 
pansion time, is given by 7 C = &irm e c/ (<jTlB 2 t), where t 
is the observer's time (neglecting the redshift correction) 
(Sari et al. 1998). After having -f m and j c , we can eas- 
ily obtain the evolution of v m and v c with time based on 
v m oc 77™ -B and v c oc 77 2 -B (see below). The remaining 
break frequency is the self-absorption one, which is given 
by v a = {5enr/[(3 - s)B~/ m ]} 3 / 5 i/ m for v a < v m < v c 
(expected at late times of the afterglow, cf. Panaitescu & 
Kumar 2000). 

The next crucial question is how the Lorentz factor 7 
decays with the observer's time because the break fre- 
quencies and the peak flux, which are needed in calcu- 
lating the observed flux, are functions of 7 and of the 
shock radius and medium density. Even if the shock is 
beamed, as long as 7 > 9$ (c s /c), the jet evolution is 
a sphere-like expansion based on the Blandford-McKee 
(1976) self-similar solution, and thus the Lorentz fac- 
tor decreases as 7 = B^E^/n*) 1 / 8 * -3 / 8 for s = 0, or 
7 = 8.8(£ 5 3M*) 1/4 £~ 1/4 for s = 2, where E 53 is the 
isotropic-equivalent energy of the jet in units of 10 53 ergs, 
and t is in units of 1 day. However, the transition of the 
jet evolution takes place at 7 ~ 9q 1 (c s /c), which in fact 
defines the break time After this time, the jet will enter 
the spreading phase with c s = c/\/3 in Rhoads (1999) or 
with c s = c in Sari et al. (1999). As a result, the Lorentz 
factor decays as 7 oc i -1 / 2 . 

After knowing the evolution of 7, we can find scaling re- 
lations of the break frequencies with time. First, we derive 
the characteristic frequency 

( £ -[3(p+2)]/[8( p -i)] 5 spherical in ISM, 
i/ m « Hp+^/Wp- 1 )], spherical in wind, (2) 

( t-(p+2)/[2( P -l)] ; jet . 



Second, the self-absorption frequency is found to evolve as 

f f [9(2- P )]/[i6(p-i)] ; spherical in ISM, 
v a <x\ t (74-4? P )/[40(p-i)] ) spherical in wind, (3) 

y t (34-19p)/[20(p-l)] ) j et _ 

Finally, the cooling frequency evolves as v c oc i -1 / 2 for a 
spherical shock in the ISM, v c oc t 1 / 2 for a spherical shock 
in the wind, and v c oc t° for a jet. 

In addition, the observed peak flux, F Vm , has been de- 
rived by many authors (e.g., Waxman 1997; Dai & Lu 
1998; Wijers & Galama 1999; Rhoads 1999; Sari et al. 
1998, 1999; Chevalier & Li 2000). In this Letter we ne- 
glect the effect of dust extinction on the peak flux because 
this effect has been discussed to be significant only for a 
highly collimated jet expanding in a dense circumstcllar 
cloud by Dai, Huang & Lu (2001). 

Therefore, we can calculate the light curves for four fre- 
quency ranges. The flux at frequencies lower than v a : 
F v<Va = Fv m {y a / 'v m ) 1/3 \v 1 V Q ) 2 , and thus evolves as 

f t (i7 P -26)/[i6(p-i)] j spherical in ISM, 
F v<Va oc I t(i3 P -i8)/[8(p-i)] j spherical in wind, (4) 

{ t[3(p-2)]/[4( P -l)] ; jet _ 

The flux above the self-absorption frequency but be- 
low the characteristic frequency is given by F Va<v<Vm = 
Fu m {v /vm) 1 / 3 , which evolves as 

C t(p+2)/[8( P -i)] j spherical in ISM, 
F Va<v<Vm ot\ f[5(2- P )]/[i2(p-i)] ! spherical in wind, 

[ t (8-6p)/[6(p-l)] } jet. 

(5) 

It is seen from equations (4) and (5) that, below i/ m , the 
light curve index is still determined by p. As a compar- 
ison, the index is independent of p in the case of p > 2 
(Sari et al. 1999). For p = 1.5 (similar to the index ob- 
tained by Malkov 1999 for Fermi acceleration in the limit 
when particles acquire a significant fraction of the shock 
energy), the flux at v < v a is approximately constant, and 
the flux at v a < v < v m increases as oc i 7 / 8 for the ISM 
case and oc t 5 ^ 12 for the wind case, respectively, as long 
as the expansion is spherical. Then, once the jet enters 
the spreading phase, the flux below the self- absorption fre- 
quency begins to decline as oc t -0 - 75 , and the flux at higher 
frequency begins to increase slowly as oc t 1 / 6 . 

If the observed high-frequency emission comes from the 
radiating electrons that are slow cooling, we have its flux 
F Vm<v<Vc = F Vm {v 1 V m )~( p_1 - 1/2 , which decays as 

C f-3(P+2)/i6 5 spherical in ISM, 
F Vm<v<Vc oc <^ t-(p+ 8 )/ 8 , spherical in wind, (6) 

{ t-(P+6)/4, jct . 

Above the cooling frequency, we obtain F V>V<1 = 
^ m (^mr (p ~ 1)/2 (i7^ c r p/2 , wh ic h declines as 

( t-(3p+io)/i6 ) spherical in ISM, 
F v>Va oc \ t-(p+6)/8 > spherical in wind, (7) 
( t-(f+6/4 7 jet. 

Bhattacharya (2001) derived light curves of the emission 
from a jet expanding in the ISM by assuming a general case 
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of 7m oc "f q in equation (1). Our light curves in equations 
(6) and (7) are consistent with his q = —1/2 result. Define 
the light curve index a and the spectral index (3 through 
F v {t) oc t~ a v~P . Table 1 summarizes the relations be- 
tween a and (3 above v m for different cases. Figures 1 and 
2 further present the a — (3 relations for 1 < p < 2 as 
well as those for p > 2, in the ISM and wind cases, respec- 
tively. We see that, for each line in these figures, the p < 2 
segment is not an extrapolation of the p > 2 segment. 

3. COMPARISON WITH THE AFTERGLOW OF GRB 010222 

The UBVRI light curve of the GRB 010222 afterglow 
has been fitted by one broken power law: F v oc t~ ai be- 
fore the break time tb and F v oc t~" 2 after if,. Here we 
summarize the light curve indices, the break time, and 
the spectral index given in the literature: (a\, ct2, tb, [3) 
are (0.60 ± 0.03, 1.31 ± 0.03, 0.48 ± 0.02 days, 1.1 ± 0.1) 
(Masettiet al. 2001), (0.80±0.05, 1.30 ±0.05, 0.72 ±0.10 
days, 0.88 ±0.10) (Stanek et al. 2001), and (0.74 ±0.05, 
1.35±0.04, 0.7±0.07days, 0.75±0.02) (Sagaretal. 2001b). 
In addition, the X-ray decay index after the break mea- 
sured by BcppoSAX is a 2 — 1.33 ± 0.04 and the spectral 
index (3 = 0.97 ± 0.05 (in 't Zand et al. 2001). A com- 
mon result of the optical and X-ray observations is that 
the light curve indeed began to steepen to ex t~ 13 about 
0.5 days after the GRB. This is the earliest observed break 
of all the studied afterglows. 

The temporal property of the afterglow from GRB 
010222 is naturally reminiscent of the jet model. Indeed, 
some authors (e.g., Stanek et al. 2001; Sagar et al. 2001b; 
Cowsik et al. 2001) attributed this afterglow to a highly 
collimated jet. Stanek et al. gave a spectral fit of their 
BVRI data, and obtained an index of (3 = 0.88 ± 0.10, 
in excellent agreement with the g'r'i'z' fit of Lee et al. 
(2001), (3 = 0.90 ± 0.03, and with the spectral index given 
by Jha et al. (2001), (3 = 0.89 ± 0.03. This implies a 
spectral index of the electron distribution, p = 2.8 in the 
slow-cooling electron regime or p = 1.8 in the fast-cooling 
electron regime. The former value of p leads to F v oc i~ 2 - 8 
at late times (Rhoads 1999; Sari et al. 1999) while the lat- 
ter value gives F v oc t~ 1,95 (see section 2). These results 
are inconsistent with the observed late-time light curve 
(oc i~ 1,3 ). Stanek et al. have noted this inconsistency. To 
save the jet model, they suggested that the spectral index 
could be intrinsically in the range of 0.5 < (3 < 0.7 due to 
the SMC-like extinction. Such a range of the spectral in- 
dex requires 2 < p < 2.4 (slow cooling) or 1 < p < 1.4 
(fast cooling). Even if the value of p becomes smaller 
for the theoretical spectral index to be compatible with 
the observed extinction-corrected spectral index, accord- 
ing to Rhoads (1999), Sari et al. (1999) and our analysis 
in section 2, we still conclude that the spreading jet model 
cannot provide an explanation for the late-time light curve 
index. We note that Sagar et al. (2001b) suggested the af- 
terglow of GRB 010222 as evidence for a highly collimated 
jet with a fast-cooling, flat-spectrum electron distribution. 
Their argument is that the emission flux from a spreading 
jet decays as oc t~ p for 1 < p < 2, which means F v oc i~ 13 
when p = 1.3, inferred by their fitting spectrum. However, 



from our analysis in section 2, we see that their argument 
is incorrect. 

An alternative explanation for the afterglow of GRB 
010222 is the expansion of a relativistic fireball or a mildly 
collimated jet in a medium with density of 10 5 — 10 6 cm~ 3 
(Masetti et al. 2001; in 't Zand et al. 2001). In 
such a dense medium, the fireball decelerated to the non- 
relativistic regime within a few days after the burst, re- 
sulting in a steepening of the light curve (Dai & Lu 1999, 
2000). in 't Zand et al. argued that the non-relativistic 
interpretation with a universal p w 2.2 value is consis- 
tent with the observations. They also noted that the 
dense-medium assumption is compatible with the observed 
rcdshift-corrected column density of <~ 2.5 x 10 22 cm~ 2 . 

4. DISCUSSION AND CONCLUSION 

We have derived light curves of the emission when a 
highly collimated jet with a flat-spectrum electron distri- 
bution (1 < p < 2) expands in the ISM or in the pre-burst 
wind. The most important finding of ours is that once 
the jet begins to spread, the light curve index becomes 
(p + 6)/4 rather than p. Therefore, the jet model appears 
to be inconsistent with the afterglow data of GRB 010222. 
ICS in the shocked medium doesn't influence the light 
curves derived in section 2. This is because for p < 2 most 
of the electron energy behind the shock front should be ra- 
diated away via both synchrotron radiation and ICS and 
thus the Compton parameter Y w (— 1 + \Jl + 4e e /es)/2 
is approximately constant (Panaitescu & Kumar 2000; Sari 
& Esin 2001). 

Two important quantities that future observations led 
by HETE-2 and Swift will provide are the light curve in- 
dex and the spectral index, which, once known, will show 
a point in Figures 1 and 2. According to the position of 
this point in these figures, one could not only obtain infor- 
mation on the dynamical evolution of a post-burst shock 
wave and the radiation regime of the accelerated electrons 
(slow cooling or fast cooling), but also infer the value of p. 

It should be emphasized that our derivations in section 
2 are based on the assumption that the electron energy 
density behind a shock is a constant fraction (e e ) of the 
total energy density of the shocked medium, as used in 
the standard afterglow shock model. If this assumption is 
invalid, the minimum Lorentz factor of the electrons, with- 
out any acceleration, could become 7 instead of equation 
(1). In such a case, the previous jet model could explain 
the afterglow of GRB 010222 if the electron energy distri- 
bution is required to be a power law with p < 2. However, 
it is unclear whether this requirement is satisfied in the 
absence of any acceleration. 

We are very grateful to the referee and D. M. Wei 
for valuable comments that significantly improved the 
manuscript, and to N. Masetti and B. Zhang for discus- 
sions. This work was supported by a RGC grant of Hong 
Kong government, the National Natural Science Foun- 
dation of China (grant 19825109), and the National 973 
Project (NKBRSF G19990754). 
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spectral index (3 
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a = 3(p + 2)/16 
a = 3(2/3 + 3)/16 


a= (p + 8)/8 
a = (2/3 + 9)/8 


(p + 6)/4 
(2/? + 7)/4 


V > v c 


[3 = p/2 


a = (3p+10)/16 
a = (3/3 + 5)/8 


a = (p + 6)/8 
a = (/3 + 3)/4 


(p + 6)/4 
(/3 + 3)/2 



Table 1 

The spectral index f3 and the light curve index a as function of p in the case of 1 < p < 2. The parameter-free 

RELATION BETWEEN a AND /3 IS GIVEN FOR EACH CASE BY ELIMINATING p. 
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Fig. 1. — A plot of the light curve index (a) versus the spectral index (]3) in the ISM case. Lines A and B correspond to a highly collimated 
but spreading jet whose (observed) high-frequency emission comes from the radiating electrons that are slow cooling {v < v c ) and fast cooling 
(y > !/ c ), respectively, and lines C and D to a spherical shock. 
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A: jet v < v c 
B: jet v > v c 
C: spherical v < v c 
D: spherical v > v 



1.0 

Spectral index (3 



2.0 



Fig. 2. — Same as Fig. 1 but in the wind case. 



